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ABSTRACT 

Recent observations suggest that some high-velocity clouds may be confined by massive dark matter halos. In 
particular, the proximity and proposed dark matter content of the Smith Cloud make it a tempting target for the 
indirect detection of dark matter annihilation. We argue that the Smith Cloud may be a better target than some 
Milky Way dwarf spheroidal satellite galaxies and use y-ray observations from the Fermi Large Area Telescope to 
search for a dark matter annihilation signal. No significant y-ray excess is found coincident with the Smith Cloud, 
and we set strong limits on the dark matter annihilation cross section assuming a spatially extended dark matter 
profile consistent with dynamical modeling of the Smith Cloud. Notably, these limits exclude the canonical thermal 
relic cross section (~3x 10 -26 cm 3 s -1 ) for dark matter masses <30 GeV annihilating via the bb or t + t _ channels 
for certain assumptions of the dark matter density profile; however, uncertainties in the dark matter content of the 
Smith Cloud may significantly weaken these constraints. 
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1. INTRODUCTION 

Many theories predict that dark matter is composed of weakly 
interacting massive particles (WIMPs), which may annihilate 
into standard model particles in regions of high dark matter 
density (e.g., Bertone et al. 2005; Baltz et al. 2008). Since 
the conventional WIMP mass scale resides in the GeV to 
TeV range, one of the principal products of dark matter 
annihilation would be a flux of y-rays. For this reason, y- ray 
searches with the Large Area Telescope on board the Fermi 
Gamma-ray Space Telescope ( Fermi-LKT ; Atwood et al. 2009) 
have provided some of the strongest limits on dark matter 
annihilation. 

Searches for dark matter annihilation traditionally target re- 
gions where an appreciable dark matter density is observation- 
ally confirmed and regions where a large dark matter density 
is strongly motivated by theory. The former category includes 
the Milky Way dwarf spheroidal satellite galaxies, where stel- 
lar velocity measurements can often constrain the dark mat- 
ter content to within a factor of two (e.g., Abdo et al. 2010; 
Geringer-Sameth & Koushiappas 2011; Ackermann et al. 
2011, 2014). The latter category includes the Galactic Center 
(e.g., Hooper & Goodenough 2011; Hooper & Linden 2011; 
Abazajian & Kaplinghat 2012; Gordon & Macias 2013; 
Gomez- Vargas et al. 2013) and clusters of galaxies (e.g., 
Ackermann et al. 2010; Sanchez-Conde et al. 2011; Ando & 
Nagai 2012; Han et al. 2012), where the predicted dark mat- 
ter signal strength is uncertain by several orders of magnitude. 
It is essential to examine multiple targets when searching for 
dark matter annihilation, since other astrophysical backgrounds 
may mimic a dark matter signal in a specific region. Addition- 
ally, the detection of dark matter in multiple Galactic environ- 
ments would provide a valuable check on the predictions of 


cosmological simulations by mapping the dark matter distribu- 
tion on Galactic and sub-Galactic scales. 

High-velocity clouds (HVCs) are a unique class of nearby 
Galactic substructures which may host significant dark matter 
content. HVCs are detected as cold clouds of neutral hydrogen 
(Hi), and are characterized by their large peculiar velocities, 
which are often incompatible with Galactic rotation (Wakker 
& van Woerden 1997). Multiple HVCs exhibit cometary mor- 
phologies, indicative of stripping by an external medium (Briins 
et al. 2000), and also have peculiar velocities which imply pre- 
vious encounters with the Milky Way disk (Mirabel & Morras 
1990). The origin of HVCs is currently unknown; however, 
several popular theories have arisen. A subset of clouds ap- 
pears to be correlated with the Magellanic Stream, suggesting 
an origin linked to interactions between the Milky Way and the 
Magellanic Clouds (Mathewson et al. 1974). Other models sug- 
gest that HVCs may be associated with globular clusters, the 
polar ring, or may be jettisoned from the Galactic nucleus. Yet 
another possibility is that HVCs may be accreted from inter- 
galactic space (Wakker & van Woerden 1997). 

The dark matter content of HVCs is controversial. In models 
where HVCs are correlated with the Magellanic Stream or are 
jettisoned from galactic environments, the dark matter content 
of HVCs is likely to be negligible. However, Blitz et al. (1999) 
argue that the spatial distribution of a sub-population of HVCs 
is consistent with the expected distribution of dark matter 
subhalos. These HVCs are expected to contain a substantial dark 
matter mass, potentially orders of magnitude larger than their H i 
mass. Braun & Burton (2000) argue that the sizes of observed H i 
clouds are not consistent with a mass profile consisting solely 
of the observed gas, and that the outer regions of the clouds 
would be unbound without a significant dark matter component. 
However, Plockinger & Hensler (2012) analyze the disruption 
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of HVCs by the Galactic disk and conclude that the majority of 
HVCs do not contain a significant dark matter component. An 
analysis of compact Hi clouds by Saul et al. (2012) find that 
the majority of HVCs are likely associated with stellar outflows 
or the Magellanic Stream, and thus possess little or no dark 
matter. Westmeier et al. (2008) note that the difference between 
the observed distribution of HVCs and the expected distribution 
of cold dark matter subhalos could be corrected if gas far from 
the Galaxy is highly ionized, preventing the detection of distant 
HVCs through 21 cm line observations. Furthermore, they note 
that only a small sub-population of HVCs must be associated 
with cold dark matter halos to match the expectations from 
simulations (Kravtsov et al. 2004). Thus, the lack of a significant 
dark matter component in some HVCs does not preclude the 
existence of a sub-population of HVCs possessing large dark 
matter content. 

The Smith Cloud is one of the better- studied HVCs, owing 
to its substantial mass, relative proximity, and location near the 
Galactic plane (Smith 1963). As a low-metallicity HVC (Hill 
et al. 2009) lacking any clear association with the Magellanic 
Stream, the Smith Cloud is a plausible candidate to host a 
significant dark matter halo — though it is worth noting that 
low-metallicity HVCs could also arise from the tidal disruption 
of metal-poor galaxies (Plockinger & Hensler 2012). The Smith 
Cloud has a striking cometary structure, indicative of a previous 
interaction with the Galactic disk. The Smith Cloud resides 
2.9 ±0.3 kpc below the Galactic plane at a Galactocentric 
distance of 7.6 zb 0.9 kpc and a heliocentric distance of 12.4 d= 
1.3 kpc. The H i content of the Smith Cloud has a projected size 
of > 3 x 1 kpc and a current mass of ~ 10 6 M© , with the brightest 
Hi emission located at l,b = 38?67, — 13741 (Lockman et al. 
2008). The three-dimensional motion of the Smith Cloud can 
be derived from its distance, morphology, and the positional 
dependence of its systemic velocity relative to the local standard 
of rest. The derived orbital parameters suggest that the Smith 
Cloud passed through the Galactic plane ~70 Myr ago and is 
predicted to cross the plane again in ~27 Myr (Lockman et al. 
2008). This observation is puzzling, since the relatively low Hi 
mass density of the Smith Cloud suggests that it should have 
been disrupted by such an encounter. Specifically, the gaseous 
component has a weak self-gravity, significantly lower than the 
ram pressure force from an interaction with the Galactic disk, 
and complete dispersion of the Smith Cloud should occur on 
timescales of less than 1 Myr. 

Recently, Nichols & Bland-Hawthorn (2009, hereafter 
NBH09) modeled the interaction of the Smith Cloud with the 
Galactic disk, employing the observed dynamical properties 
of the Smith Cloud as determined by Lockman et al. (2008). 
By modeling the tidal disruption and ram-pressure stripping of 
the Smith Cloud during its collision with the Galactic plane, 
NBH09 determined that an additional massive component (i.e., 
dark matter) is required for the survival of the gas cloud. To 
address uncertainties in the dark matter density profile, NBH09 
simulated situations where the dark matter density follows a 
Navarro-Frenk-White (NFW; Navarro et al. 1996), Burkert 
(Burkert 1995), or Einasto (Springel et al. 2008) profile. To 
constrain uncertainties resulting from prior encounters of the 
Smith Cloud with the Galactic plane, NBH09 considered sce- 
narios in which the Smith Cloud has transited the Galactic plane 
only once, as well as scenarios where the gas distribution of the 
Smith Cloud has reached a steady state due to repeated dynam- 
ical interactions with the plane. In all cases, NBH09 found that 
a significant dark matter component is required to prevent the 


disruption of the Smith Cloud. To gravitationally bind the Smith 
Cloud, NBH09 calculated that the dark matter mass within 1 kpc 
of the cloud center must exceed 2 x 10 8 M© prior to the most 
recent stripping event. Moreover, tidal stripping during the most 
recent interaction with the Galactic plane is expected to have 
decreased the dark matter mass by approximately a factor of 
two, compared to a factor of five decrease in the gas content. 
This indicates that the Smith Cloud is currently composed of 
nearly 99% dark matter, comparable to the dark matter fraction 
in dwarf spheroidal galaxies. 

In this work, we use Fermi-LAI data to search for excess 
y- ray emission coincident with the Smith Cloud, accounting 
for the y-ray emission from nearby point-like sources, Galac- 
tic foregrounds, and extragalactic backgrounds. Conventional 
models of the Galactic y-ray foreground include a component 
representing the emission from Galactic cosmic rays interacting 
with interstellar gas in the Smith Cloud. It is essential to remove 
this component from the foreground model before searching for 
residual y-ray emission associated with the Smith Cloud. No 
substantial residual flux is observed coincident with the Smith 
Cloud, and we place strong constraints on the dark matter an- 
nihilation cross section in scenarios where the Smith Cloud is 
dynamically bound by a dark matter halo with the character- 
istics derived by NBH09. Interestingly, we note that in many 
regimes these limits are stronger than those obtained from ob- 
servations of dwarf spheroidal galaxies due to the proximity of 
the Smith Cloud. However, we stress that the dark matter con- 
tent of the Smith Cloud is far more uncertain than that of dwarf 
spheroidal galaxies. 

The outline of this paper is as follows. In Section 2 we model 
the y-ray emission from dark matter annihilation in the Smith 
Cloud. In Section 3 we summarize the procedure for creating 
a Galactic foreground model that excludes the Smith Cloud. 
In Section 4 we describe our analysis of the Fermi-LAI data 
in the region of the Smith Cloud. In Section 5 we set limits 
on the y-ray flux associated with the Smith Cloud and constrain 
the dark matter annihilation cross section. Finally, in Section 6 
we discuss the dependence of these limits on the dark matter 
content of the Smith Cloud, compare these limits with those 
derived from observations of dwarf spheroidal galaxies and the 
Galactic Center, and discuss how this analysis could be extended 
to a larger population of HVCs. 


2. DARK MATTER MODELS OF THE SMITH CLOUD 


The integrated y-ray flux at Earth, cj) s ( ph cm -2 s -1 ), expected 
from dark matter annihilation in a density distribution, p(r), is 
given by 
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Here, the <D PP term depends on the particle physics properties 
of dark matter — i.e., the thermally averaged annihilation cross 
section, (crv), the particle mass, m P M, and the differential 
y-ray yield per annihilation, dN y /dE y , integrated over the 
experimental energy range from E m i n to £ max . The J-factor is 
the line-of- sight integral through the dark matter distribution 
integrated over a solid angle, AO. Qualitatively, the J-factor 
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Table 1 

Summary of Smith Cloud Dark Matter Halo Parameters 


Profile 

r s 

(kpc) 

PO 

( Mq kpc -3 ) 

A/tidal 

(Mq) 

J-factor 

( GeV 2 cm -5 sr) 

NFW 

1.04 

3.7 x 10 7 

1.1 x 10 8 

9.6 x 10 19 

Burkert 

1.04 

3.7 x 10 7 

1.3 x 10 8 

4.2 x 10 18 

Einasto 

1.04 

9.2 x 10 6 

2.0 x 10 8 

1.8 x 10 2 ° 


Note. Integrated J-factors are calculated over a solid-angle cone with radius 1° 
(AQ ~ 9.6 x 10 -4 sr). 

encapsulates the spatial distribution of the dark matter signal, 
while Opp sets its spectral character. 

There is significant uncertainty in the dark matter density 
profile of the Smith Cloud, thus we calculate the J-factor for 
each of the three dark matter halo profiles fit by NBH09. The 
parameters of each profile were derived from the total dark 
matter tidal mass, M^i, required to confine gas in the Smith 
Cloud during its most recent interaction with the Galactic disk. 
These tidal masses are calculated independently by NBH09 for 
each dark matter profile. We describe each of the dark matter 
density profiles in terms of a scale radius, r s , and a scale density, 
Po, as listed in Table 1. The Einasto profile depends on an 
additional parameter a which is set to a value of 0.17: 


p(r) = 

P0 r s 
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To avoid peripheral regions where tidal stripping may alter the 
dark matter density, we truncate our model of the y-ray intensity 
profile 1° from the center of the Smith Cloud. To simplify 
comparisons with other dark matter annihilation targets (i.e., 
dwarf spheroidal galaxies), we compute the integrated J-factor 
from the Smith Cloud within this 1° radius (Table 1). This radius 
contains ~60% of the total predicted y- ray flux when cuspy 
NFW or Einasto profiles are assumed and ~10% of the total 
predicted flux from the cored Burkert model. Thus, this choice 
of radius yields a conservative estimate for the total J-factor of 
the Smith Cloud since the dark matter distribution may extend 
to several degrees. 

3. GALACTIC FOREGROUND MODELING 

The observed foreground y-ray emission from the region 
surrounding the Smith Cloud is dominated by 7r° -decay emis- 
sion produced from cosmic rays interacting with the atomic 
and molecular hydrogen gas of the Milky Way. 8 The GALPROP 
cosmic-ray propagation code can be used to model the dif- 
fuse Galactic y-ray foreground from processes such as inelastic 
hadronic collisions, bremsstrahlung, and inverse-Compton scat- 
tering. 9 * GALPROP accounts for effects such as diffusion, reac- 
celeration, and energy loss via mechanisms such as synchrotron 
radiation (Strong & Moskalenko 1998; Strong et al. 2009). 


8 The y - ray emission from inelastic hadronic interactions is composed of 
many processes, the most important of which being the production of tt°, 
which decays primarily to yy. 

9 http://galprop.stanford.edu 


The distribution of target material is derived from surveys of 
the 2.6 mm CO and 21 cm Hi lines, supplemented with inter- 
stellar reddening maps from infrared observations of interstellar 
dust. Notably, the official Fermi-LXY model of Galactic diffuse 
emission recommended for discrete source analysis includes a 
y-ray emission component associated with the Smith Cloud. 10 
We remove gas correlated with the Smith Cloud from our 
analysis for two reasons. First, the intensity and spectrum of 
cosmic rays are poorly constrained at the distance of the Smith 
Cloud, which leads to considerable uncertainty in the predicted 
y- ray flux. Second, removing gas from the Smith Cloud elimi- 
nates a potentially degenerate emission component which may 
result in artificially strong limits on the dark matter annihilation 
rate within the cloud. 

We create Galactocentric annuli for the H i gas distribution 
by transforming 21 cm brightness temperatures into column 
densities using the composite LAB survey (Kalberla et al. 2005) 
and the Galactic rotation curve given by Clemens (1985). We 
follow the procedure employed by Ackermann et al. (2012) to 
excise the gas associated with the Smith Cloud from the Galactic 
gas distribution. Specifically, we remove gas in the region 
36° ^ l ^ 46° and —20° ^ b ^ — 10°, which has a velocity with 
respect to the local standard of rest in the range 70-125 km s _1 
(Figure l). 11 The primary uncertainty in the conversion from 
brightness temperature to column density comes from the 
assumed spin temperature (T s ) used to correct for the opacity 
of the 21 cm line. We find that the gas density in the region 
of the Smith Cloud changes by <15% when the assumed spin 
temperature is changed from T s = 125 K to Ts = 10 5 K (i.e., 
the gas is optically thin). When analyzing the y- ray data we 
set T s = 125 K and find that this choice has little impact 
on our results. We follow the procedure of Ackermann et al. 
(2012) to trace the CO distribution from the 2 mm composite 
survey of Dame et al. (2001). Due to the small CO content of 
HVCs (Akeson & Blitz 1999), we do not alter the CO map 
in the region of the Smith Cloud. The Galactic foreground 
also contains a contribution from neutral gas that cannot be 


10 http://fermi.gsfc.nasa.gov/ssc/data/access/lat/BackgroundModels.html 

11 Gas with velocity > 125 km s -1 contributes less than 0.8% of the total 
column density. 



-100 -50 0 50 100 

Velocity (km s -1 ) 

Figure 1. Latitude-velocity distribution of Galactic H i gas from the LAB survey 
(Kalberla et al. 2005) integrated over the longitude range of the Smith Cloud 
(36° ^46°). The color represents the integrated brightness temperature of 

the 21 cm H i line as a function of latitude and velocity with respect to the local 
standard of rest. Gas associated with the Smith Cloud is enclosed by the black 
box and is removed from our Galactic foreground model. 

(A color version of this figure is available in the online journal.) 
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Figure 2. 15° x 15° ROI surrounding the Smith Cloud in the energy range from 500 MeV to 500 GeV. The gray contours represent the Hi column density associated 
with the Smith Cloud (1 x 10 20 cm -2 < Nhi < 2.7 x 10 20 cm -2 ), while the over-plotted circle shows the 1° truncation radius for the dark matter profile. Left: 
observed y- ray counts map smoothed by a Gaussian kernel with standard deviation 0?1. Center: map of the background y-ray emission model fit to the Fermi -LAI 
data including diffuse and point-like backgrounds. Right: the Poisson probability of finding the observed number of counts in each pixel given the model prediction 
expressed as a Gaussian significance. 

(A color version of this figure is available in the online journal.) 



traced by the combination of Hi and CO (so-called dark gas). 
We follow the procedure of Ackermann et al. (2012) to trace 
the dark gas using the E(B — V) reddening maps of Schlegel 
et al. (1998). We incorporate a dark gas correction into the 
H i map after the Smith Cloud has been removed (Ackermann 
et al. 2012). 12 We note that our procedure for removing the gas 
content of the Smith Cloud is very similar to the method used to 
remove gas associated with the Magellanic Clouds and M3 1 (see 
Appendix B of Ackermann et al. 2012). 

These observations of the Milky Way gas profile supple- 
mented by infrared observations of Galactic dust are input 
into the GALPR0P code to model the diffuse y- ray emission 
corresponding to hadronic collisions, inverse-Compton scat- 
tering and bremsstrahlung radiation. To provide an accurate 
model for diffuse emission in the region of the Smith Cloud, 
we adopt the best-fit propagation parameters given by Trotta 
et al. (2011), specifically a convectionless diffusion constant of 
8.32 x 10 28 cm 2 s -1 at a momentum of 4 GeV, with a power- 
law momentum scaling D(p) a p 0,3 *, a height for the diffusion 
region of 5.4 kpc, and an Alfven velocity of 38.4 km s _1 . These 
parameters were inferred from a Bayesian analysis including the 
Fermi-LAI data as an input, and the resulting model is well-fit 
to the Galactic diffuse y-ray emission at intermediate latitudes 
corresponding to the Smith Cloud. From this model we produce 
energy-dependent maps for the y- ray emission from hadronic 
emission, bremsstrahlung, and inverse-Compton scattering. In 
principle, we would consider any alterations to the propaga- 
tion parameters which are consistent with the local cosmic-ray 
primary-to-secondary ratios measured by satellite and balloon 
experiments. However, we find that this first attempt yields an 
accurate model of the observed diffuse y-ray emission in the 
region of the Smith Cloud and no additional parameter- space 
scan is necessary. 

4. DATA ANALYSIS 

To search for excess y-ray emission coincident with the 
Smith Cloud, we select a data sample corresponding to the 
first five years of Fermi-LAI operation (2008 August 4 to 2013 


12 The E(B — V) correction excludes the Smith Cloud due to its low 
metallicity. 


August 4). We select events from the P7REP CLEAN class in 
the energy range from 500 MeV to 500 GeV and within a 15° 
radius of the Smith Cloud (/, b = 3 8? 67, — 13241). Extending 
this analysis to lower energies would translate to a minor 
improvement in the sensitivity to low-mass dark matter models; 
however, below 500 MeV the rapidly changing effective area 
results in a stronger dependence on the spectral model assumed 
for the Smith Cloud. To reduce y-ray contamination from the 
Earth’s limb, we reject events with zenith angles larger than 100° 
and events collected during time periods when the magnitude of 
the rocking angle of the Fermi-LAI was greater than 52°. 

We perform a binned maximum likelihood analysis of the 
15° x 15° region-of-interest (ROI) surrounding the Smith Cloud 
(Figure 2). We bin the Fermi-LAI data in this ROI into 0? 1 pixels 
and 24 logarithmically spaced bins of energy from 500 MeV to 
500 GeV. We model the diffuse emission in this region using 
the templates for the hadronic, bremsstrahlung, and inverse- 
Compton emission derived in the previous section. Because the 
hadronic and bremsstrahlung components are morphologically 
similar (both trace the interstellar gas in the Milky Way), we 
tie their relative normalizations in the y-ray fit. In addition to 
the diffuse Galactic foregrounds, the y-ray data includes an 
isotropic contribution from extragalactic y-rays and charged 
particle contamination. The spectrum of the isotropic y-ray 
background is routinely derived from a high-latitude (\b\ > 10°) 
fit to the Fermi-LAI data, and is therefore dependent on the 
data selection and on the modeling of other y-ray emission 
components (i.e., the Galactic diffuse emission). It is difficult to 
derive the detailed spectrum of this component locally in the ROI 
of the Smith Cloud due to limited statistics at high energies and a 
morphological degeneracy with the inverse-Compton emission. 
Thus, we model the spectrum of the isotropic component with 
a broken power-law model which is simultaneously fit to the 
Fermi-LAI data in the Smith Cloud ROI. While a broken power- 
law model offers a reasonable fit to the Smith Cloud ROI, it 
does not capture the detailed energy dependence of the residual 
background. To quantify the impact of this simplification we 
also perform the analysis with the standard isotropic background 
model, iso_clean_v05.txt , 13 and find that the results change by 


13 http://fermi.gsfc.nasa.gov/ssc/data/access/lat/BackgroundModels.html 
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Figure 3. Bin-by-bin energy-flux upper limits and expected sensitivities at 95% 
CL for the Smith Cloud assuming an NFW dark matter profile. The 95% CL 
upper limits derived from the data are shown by the black arrows. The median 
sensitivity is shown by the dashed black line while the 68% and 95% containment 
regions are indicated by the shaded bands. 

(A color version of this figure is available in the online journal.) 

<15%, which is much smaller than the fractional uncertainty 
in the dark matter distribution of the Smith Cloud. In addition 
to the diffuse components, our model includes all sources from 
the second Fermi-LAI source catalog within 20° of the Smith 
Cloud (Nolan et al. 2012). The flux normalizations of sources 
within 6° of the Smith Cloud are left free in the fit. 

Since we are specifically interested in y-ray emission as- 
sociated with dark matter annihilation in the Smith Cloud, we 
model the Smith Cloud itself according to the spatially extended 
dark matter profiles described in Section 2. We follow the pre- 
scription of Ackermann et al. (2014) to perform a bin-by-bin 
likelihood analysis of the y-ray emission coincident with the 
Smith Cloud. We first perform a single fit over the entire energy 
range to fix the normalizations of the diffuse and point- like back- 
ground sources. 14 We then perform a likelihood scan over the 
normalization of the putative Smith Cloud dark matter source in- 
dependently in each energy bin (this procedure is similar to that 
used to evaluate the spectral energy distribution of a source). By 
analyzing each energy bin separately, we avoid selecting a sin- 
gle spectral shape to span the entire energy range at the expense 
of introducing additional parameters into the fit. The common 
spectral model-dependent likelihood can be reconstructed by 
tying the signal normalization across the energy bins. As a con- 
sequence, computing a single bin-by-bin likelihood function 
allows us to subsequently test many spectral models rapidly. 
The bin-by-bin likelihood is additionally powerful because it 
presents the y- ray data in a way that makes minimal assump- 
tions about the y- ray spectrum of the source of interest. While 
the bin-by-bin likelihood function is essentially independent of 
spectral assumptions, it does depend on the spatial model of the 
Smith Cloud and must be derived for each profile in Table 1 . 

5. RESULTS 

We find no statistically significant excess y- ray emission co- 
incident with the Smith Cloud in the energy range from 500 MeV 
to 500 GeV, and we set 95% confidence level (CL) upper limits 
on the y- ray flux within each energy bin (Ligure 3). These lim- 
its agree well with the expected sensitivity as calculated from 


14 Fixing the normalizations of the background sources at their globally fit 
values avoids poor convergence in the fitting procedure resulting from the fine 
binning in energy and the degeneracy of the diffuse background components at 
high latitude. 



Figure 4. Upper limits at 95% CL on the dark matter annihilation cross section 
as a function of the dark matter particle mass assuming the best-fit NFW 
profile from Table 1 and annihilation through the bb, x + x ~ , /z + /z - , and W + W~ 
channels. 

(A color version of this figure is available in the online journal.) 



mvM (GeV) 

Figure 5. Upper limits at 95% CL for the dark matter annihilation cross section 
to bb as a function of the dark matter particle mass for different assumptions 
of the dark matter density profile. Each profile is normalized to the best-fit 
dark matter halo mass calculated by NBH09 (Table 1). The annihilation signal 
is truncated at a radius of 1° from the assumed center of the Smith Cloud to 
mitigate possible impacts from tidal stripping. 

(A color version of this figure is available in the online journal.) 

500 simulations of Fermi-LAI instrument performance and the 
background y-ray sources in the Smith Cloud ROI. Under the 
assumption that the Smith Cloud is confined by a halo of dark 
matter, we utilize the absence of a y-ray signal to set constraints 
on the dark matter annihilation cross section. Assuming an NEW 
profile with parameters listed in Table 1 , we place constraints on 
the cross section for dark matter particles annihilating through 
the bb, t + t~, /x + /x _ , and W + W~ channels (Ligure 4). 15 As- 
suming an NEW profile, these constraints exclude the canonical 
thermal relic cross section (~3 x 10 -26 cnr 5 s _1 ) for dark mat- 
ter masses <30 GeV annihilating via the bb or t + t~ channels. 
However, the J-factor calculated for the inner 1° of the Smith 
Cloud varies by more than an order of magnitude depending on 
the shape of the assumed dark matter profile. Current observa- 
tions of the Smith Cloud offer no constraints on the shape of 
its dark matter profile and this uncertainty propagates directly 
into the constraints on the dark matter annihilation cross section 
(Ligure 5). Uncertainty in the shape and content of the Smith 


15 Dark matter annihilation spectra were generated using DMFIT (Jeltema & 
Profumo 2008; Ackermann et al. 2014). 
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Cloud dark matter halo dominates statistical uncertainties in 
the y- ray data, systematic uncertainties in the modeling of the 
Fermi-LAI instrument response, and systematic effects arising 
from the modeling of the diffuse y-ray backgrounds. 

The most significant deviation from the background- only 
hypothesis arises when fitting a dark matter particle with 
mass 5 GeV annihilating to z + z~. Incorporating this additional 
component increases the log-likelihood slightly, corresponding 
to a test statistic (TS) of TS = 4.7. This deviation is well 
below the conventional threshold (TS > 25) for the detection 
of discrete y - ray sources (Nolan et al. 2012) and is spectrally 
consistent with hadronic emission produced from cosmic-ray 
interactions with the gas of the Smith Cloud. This hadronic 
component was explicitly removed from the astrophysical 
background model to provide a conservative constraint on the 
dark matter annihilation cross section. If the gas associated 
with the Smith Cloud is not removed from the astrophysical 
background model, the significance of the additional dark matter 
component is reduced to TS ~ 2. The lack of significant 
cosmic-ray induced y-ray emission from the Smith Cloud is not 
surprising due to the mass and distance of the cloud; however, 
a more detailed investigation is beyond the scope of this work. 

6. DISCUSSION AND CONCLUSIONS 

Our analysis found no statistically significant y-ray excess 
coincident with the Smith Cloud, allowing us to place limits on 
the dark matter annihilation cross section. Assuming the best- 
fit NFW profile from NBH09, these constraints are stronger 
than those derived from a four-year analysis of 15 stacked 
dwarf spheroidal galaxies for dark matter masses <lTeV 
(Ackermann et al. 2014). The strength of these limits stems 
primarily from the fact that the Smith Cloud has a nominal 
J-factor approximately a factor of two higher than any of 
the dwarf spheroidal galaxies. Thus, we conclude that HVCs 
may be excellent indirect detection targets, motivating further 
investigation of the dynamics of these systems. 

It is important to note that while the constraints from dwarf 
spheroidal galaxies incorporate the uncertainty in the dark mat- 
ter content of these objects, a similarly detailed understanding 
of the Smith Cloud is currently lacking. In Figure 5, we show 
that the limits on the dark matter annihilation cross section cal- 
culated under the assumption of a Burkert density profile are 
weaker than those calculated for the NFW and Einasto profiles 
by a factor of ~40. It is likely that none of these analytic profiles 
realistically describes the dark matter profile over the full ex- 
tent of the Smith Cloud, and thus an accurate description of the 
Smith Cloud density profile (including its overall normalization) 
currently remains the largest uncertainty in setting constraints 
on the dark matter annihilation cross section. 

We emphasize that identifying additional targets for indirect 
detection is particularly important in light of recent reports 
of excess y- ray emission compared to current y- ray diffuse 
emission models from the direction of the Galactic Center, 
which are possibly consistent with dark matter annihilation 
(e.g., Gordon & Macias 2013; Abazajian et al. 2014; Daylan 
et al. 2014). The cross section fit to the proposed Galactic 
Center excess barely evades the constraints derived from the 
Smith Cloud assuming an Einasto dark matter profile but lies 
a factor of ~40 below the limits derived assuming a Burkert 
profile. We find no statistically significant excess of y-rays in 
the 1-3 GeV region where the y- ray excess from the Galactic 
Center is most pronounced. Even if a statistically significant 
y- ray excess were to be found coincident with the Smith Cloud, 


more thorough analyses of diffuse astrophysical y- ray emission 
would be necessary in order to determine the origin of the 
emission. We emphasize that new indirect detection targets 
offer a compelling and complementary method to improve 
our sensitivity to dark matter annihilation using currently 
available data. 

The techniques developed in this paper can be readily ex- 
tended to other Milky Way HVCs. More than 560 HVCs have 
been detected surrounding the Milky Way (Wakker & van 
Woerden 1991). While the distances to many HVCs are un- 
certain, these systems are relatively close to the solar position 
compared to the population of dwarf spheroidal galaxies. Fu- 
ture efforts to understand the dark matter content of HVCs and 
subsequently to perform a joint likelihood analysis of multiple 
HVCs may provide a sensitive test of dark matter annihilation. 

We note that during the final preparation of this paper, Nichols 
et al. (2014) reported upper limits on the dark matter annihilation 
cross section from Fermi-LAI observations of the Smith Cloud. 
While our results are qualitatively similar, we note three key 
differences between the Fermi-LAI analyses in these studies. 
First, we remove gas correlated with the Smith Cloud before 
producing constraints on the Smith Cloud y- ray emission, thus 
setting conservative constraints on the dark matter annihilation 
signal. Second, we utilize detailed spectral information when 
calculating upper limits on the dark matter annihilation cross 
section, rather than an integrated flux upper limit derived 
assuming a fixed power-law model for the y- ray emission 
from the Smith Cloud. This allows us to more sensitively test 
specific dark matter masses and annihilation channels. Third, 
we model the dark matter content of the Smith Cloud as a 
spatially extended source with a distribution that extends out 
to 1°, consistent with the profiles reported in NBH09. This 
last distinction significantly enhances the predicted dark matter 
annihilation signal from the Smith Cloud, and improves the 
constraints on the dark matter annihilation cross section by more 
than an order of magnitude (for an NFW profile) compared to 
the analysis of Nichols et al. (2014). 
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